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We show that the time evolution of entanglement under incoherent environment coupling can be 
faithfully recovered by monitoring the system according to a suitable measurement scheme. 
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Quantum entanglement, arguably the most non- 
classical feature of quantum mechanics, and considered 
as one of the key resources for quantum information 
processing, remains a puzzle for our intuition. Gen- 
eral characterizations of the static entanglement and/or 
the separability properties of a given quantum state re- 
fer to the intricate geometry of tensor spaces, or to the 
operation of linear operations thereon. No general ob- 
servable is known so far which would complement such 
essentially mathematical notions with a specific exper- 
imental measurement setup. Here we come up with a 
radically different, dynamical characterization of entan- 
glement, through the continuous observation of a quan- 
tum system which evolves under incoherent coupling to 
an environment. We show that there is an optimal mea- 
surement strategy which provides an optimal decompo- 
sition of the time evolved, mixed system state, in the 
following sense: The entanglement of the time evolved 
mixed state is given by the average over the pure state 
entanglement of the single realisations of the stochastic 
time evolution. The weight of the different pure states is 
herein determined by their relative abundance during the 
detection process. Since arbitrary measurement prescrip- 
tions in general do not generate mixed state decomposi- 
tions which are optimal in this sense, our finding implies 
that there are specific system observables which filter out 
the minimal amount of nonclassical correlations inscribed 
into a dynamically evolved, mixed quantum state. 

Consider a bipartite quantum system composed of sub- 
systems A and -B, interacting with its environment. Due 
to this coupling, an initial pure state of the compos- 
ite system will evolve into a mixed state in a way 
governed by the master equation p = LkP- The su- 
peroperators Ck describe the effects of the environment 
on the system, and usually take the form 



decomposition of p, 



JlJkp- pJlJk 



(1) 



where the operators Jk depend on the specific physical 
situation under study. 

To extract the time evolution of entanglement under 
this incoherent dynamics, one solution is to evaluate a 
given entanglement measure M{p) for the solution p{t), 
at all times t. One starts from one of the known pure 
state measures M(^') together with a pure state 



(2) 



where the pi are the positive, normalized weights of each 
pure state The most naive generalization for a 

mixed state would then be to consider the average 



(3) 



which, however, is not uniquely defined, since the de- 
composition 121 is not unique. Therefore, the proper def- 
inition of M{p) is the infimum of all possible averages 
M 4]. Finding this infimum has two main drawbacks: 
(i) it turns into a hard numerical problem for higher di- 
mensional or multipartite systems, and, (ii) even for bi- 
partite qubits, where analytical solutions for some mea- 
sures M{p) are known, there is no obvious general inter- 
pretation of the optimal decomposition, which gives the 
infimum, in physical terms. 

Our approach here will be to avoid the direct use of 
the mixed state solutions by substituting them by 
physically motivated ensembles of pure states. To do 
this, instead of solving the evolution equation for the den- 
sity operator, we will follow a stochastic time evolution 
d H 13 of the initially pure state. This combines ran- 
domly occurring quantum jumps, defined by the action 
of the operators Jk in Eq. ^ , with periods of continuous 
evolution, generated by a non-hermitian effective Hamil- 
tonian. The resulting quantum trajectories, see Fig. ^ 
are known to provide, after averaging over many inde- 
pendent realisations, the same result as the master equa- 
tion. Moreover, a single trajectory can be understood 
as the experimentally accessible information through the 
continuous monitoring of the system |E H, 13 ■ In other 
words, one possible decomposition of the form ^ is given 
by the set of the possible outcomes of different runs 
of one and the same experiment, after a monitoring time 
t. 

Note, however, that the stochastic time evolution in- 
duced by the jump operators Jk does not necessarily yield 
the optimal decomposition which minimises M. Yet, 
also the choice of the jump operators is not unique: 
Indeed, the definition of new jumps as linear combina- 
tions of the original ones, Ju = UkiJi, with U a uni- 
tary matrix, or by the addition of a complex number /z. 
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FIG. 1: Simplified schematic representation of stochastic time 
evolutions generated by continuous monitoring of the system 
dynamics through a jump operator J: During each time step 
of (short) duration 5t, the system evolves with a certain prob- 
ability 5p into the state l^'Jt), and with probability (1 — 5p) 
into the state I'I'it), depending on the result of the action of 
J on the initial state |^'o)- Upon iterative application of this 
method, a large number of different states can be accessed, 
giving rise to different stochastic pure state trajectories. At 
a given time t — nSt, this process generates a pure state de- 
composition of the time evolved density matrix p{t). 

J'j^ = (n ± Jk) I \/2, leaves Eq. |^ unaltered. More gen- 
erally, one can perform both changes simultaneously, to 

obtain the general jump operators £fc,± = ^/i ± Jfe^ /\/2. 

Since each such choice of jump operators produces a 
different unravelling of the same master equation (J^l, i-e., 
a different decomposition of p{t), we have to search for 
that type of monitoring which reproduces the desired op- 
timal decomposition, upon averaging over the different 
realisations at time t. We claim that such a measure- 
ment prescription can always be found, for a given initial 
state, and given environment coupling. While a rigor- 
ous mathematical proof of this conjecture is still missing, 
we now present strong evidence in its favour, for various 
decoherence models. 

We focus on the dynamics of bipartite two-levels sys- 
tems, where the results obtained by our method can read- 
ily be compared with known analytical solutions 0, 0| . 
We assume that each subsystem interacts independently 
with its own environment, and consider those situations 
where decoherence is induced by dissipation (zero tem- 
perature reservoir), by dephasing, or by noise (infinite 
temperature reservoir). For simplicity, we set all de- 
cay rates Ffc = F in Eq. We start with the gen- 
eral jump operators Lk^± as defined above, and search 
for the ones which yield the optimal solution. The first 
step in such procedure is to follow all possible states 
|\l/5t) of the system, at a short time St after the ini- 
tial time ^0 = 0; and to calculate the average mea- 
sure M{St). Since the possible states l^'^j) which can 
be reached after this short time interval are given by 
the action of the N jump operators, plus the state 
l^5t^^) which results from continuous evolution (see fig- 
ure^, the expression for the average measure is simply 
Mm = (1 - Eti Spk) M(vI/^+i) + SpkMi^i), 
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FIG. 2: Entanglement dynamics of a bipartite quantum state 
initially prepared in the maximally entangled state |v|/o) = 
(1 00) -|- |ll))/\/2, under incoherent coupling to a zero temper- 
ature environment. AI is here derived from the entanglement 
of formation .2]. The symbols result from an average over 
different realisations of the stochastic time evolution under 
the action of the general jump operators Lk,±, while the solid 
line represents the exact result derived from the integration of 
the master equation. Open squares {fi = 1) and filled circles 
{jj, = 3) are obtained by progressively improved unravellings, 
with an optimal choice of the unitary U. Open triangles in- 
dicate the result of an arbitrarily chosen unravelling, with no 
optimization, neither of U nor of p.. 

and can be minimised over the different possible unravel- 
lings. This ensures that the short time behaviour of the 
average measure is optimal. 

In fact, for many initial states, such as the Bell states 
(|00) + |ll))/\/2 and (|01) + |10))/\/2 under dephasing, 
or a state |*o) = "^oolOO) -I- -(/"oilOl) -I- -^lollO) coupled 
to a zero temperature environment, the minimisation for 
the initial step can be shown to be not only a necessary 
but also a sufficient condition, for M{t) to coincide with 
M{p{t)). This is a striking feature, since the optimisa- 
tion needs to be performed only once, for the initial time 
step, but leads to a faithful, optimal time evolution for 
arbitrary t > t^. 

In other cases, however, the correct behaviour at short 
times is not sufficient, and further search for the op- 
timal unravelling is needed, though now in a reduced 
parameter space (spanned by the parameters of ifc,±). 
Such search is illustrated in Fig. |21 for a Bell initial state 
|4'o) = (|00) + |ll))/\/2, coupled to a zero temperature 
reservoir, and for M{^!) given by the familiar entangle- 
ment of formation 2\. The condition for the correct ini- 
tial decay fixes most of the parameters of the unitary ma- 
trix J7, and enforces /i > 1/^/2. The open squares and 
filled circles correspond to unravellings satisfying these 
conditions for /i = 1, and = 3, respectively, while tri- 
angles are obtained by a non-optimal choice of U and /x. 
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FIG. 3: Time evolution of the stochastically generated en- 
tanglement of formation (which we use here for the explicit 
evaluation of M), during the application of a CNOT gate 
in a two qubit system coupled to a dephasing environment. 
Adequate choice of the jump operators L^^-i- leads to perfect 
agreement (circles) with the exact result (solid line) during 
the gate time. A non-optimal choice of unravelling leads to 
large deviations from the exact solution. 



Although both values of /i induce the correct behaviour 
for sufRciently short times, only the larger value of /j, 
gives good agreement for long times. The inset shows 
the situation for the same initial state, though coupled 
to an infinite temperature environment. The unravelling 
is chosen to give the correct initial slope, but the agree- 
ment is also very good for the asymptotics. The only dis- 
crepancy persists at intermediate times, where the exact 
solution for M{p(t)) exhibits a discontinuous derivative. 

We can allow for more general dynamics, with a non- 
trivial unitary evolution in the presence of incoherent en- 
vironment coupling. This is relevant when it comes to 
describe, for example, gate operations performed in open 
quantum systems. To illustrate such a situation, we sim- 
ulated the action of a CNOT gate, in a system coupled 
to a dephasing reservoir as shown in Fig. Starting 
with a separable initial state |^'o) = (|00) -I- |10))/V2, 
the gate would generate the maximally entangled state 
= (|00) + |11))/V2, after a time Tgato, chosen to be 
five times faster than the decoherence time scale 1 /F (see 
Eq. Q). However, due to dephasing, the entanglement 
at Tgatc will not be maximal, since the target state prepa- 
ration will not be perfect. As in our above examples one 
can find an unravelling which generates an optimal pure 
state decomposition, in excellent agreement with the ex- 
act entanglement evolution obtained from the solution 
p{t). Once again, a non-optimal unravelling suggests a 
very different entanglement dynamics, with apparently 
almost perfect gate performance - in pronounced con- 
trast to the correct result. 

Finally, let us consider the example of a tripartite qubit 



FIG. 4: Multipartite entanglement decay for a GHZ state un- 
der dephasing dynamics. Filled circles represent the tripar- 
tite concurrence C3(t) (substituted in the explicit evaluation 
of M) Q obtained after averaging over 1000 realisations of 
the unravelling defined by the jump operators — 
[k = 1, 2, 3) (which actually coincide with the coupling oper- 
ators in the dephasing master equation), in very good agree- 
ment with the exact solution indicated by the solid line. A 
non-optimal unravelling (open circles) considerably overesti- 
mates the actual mixed state entanglement of the time evolved 
state. 



system, initially prepared in a GHZ or in a W state, and 
in contact with dephasing or zero temperature environ- 
ments, respectively. In both cases, the optimal unravel- 
ling produces excellent agreement with the exact result 
for the tripartite concurrence C3 it) 0, |n| , as illustrated 
in Fig. 21 for a GHZ initial state subject to dephasing. 
In this exemplary case, the optimal detection is pro- 
vided by the original jump operators Jk = cr^'' a^^^ , with 
k G {1,2,3}. Again, a non-optimal unravelling leads to 
very different predictions. 

To conclude, we have shown that mixed state entan- 
glement can be unraveled by a suitable measurement pre- 
scription. The latter can be derived from the initial state 
of the composite quantum system under study, together 
with the specific type of environment coupli ng. There is 
thus no ambiguity (which was suggested in |l2|) in the 
definition of mixed state entanglement with respect to 
different unravellings: the minimal nonclassical correla- 
tions needed to characterize a dynamically evolving quan- 
tum state are idenitified by an optimal (in general not 
unique) monitoring prescription, and any in this sense 
non-optimal unravelling provides a non-optimal estimate 
thereof. 

Surprisingly, in all examples considered, we so far al- 
ways succeeded to come up with a time-independent op- 
timal unravelling, i.e., the monitoring strategy does not 
need to be readapted during the incoherent time evolu- 
tion. This suggests a kind of uniform convergence of the 
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stochastic average towards the exact sokition (which fol- 
lows from the evaluation of the convex roof of a given 
entanglement measure on the time evolved density ma- 
trix), for the optimal unravelling, and implies a subtle 
relation between the temporal evolution of entanglement 
and the ensemble of pure states which can be reached 



during the single realisations generated by the optimal 

jump operators. 
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